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SUMMARY 

The nature of nonadiabatio, friotionless , steady flow of a com- 
pressible fluid. In a pipe of constant cross section is summarized.. 
The flow conditions can most conveniently he described in terms of 
the local Mach number M. It was found that steady flow cannot 
oocur when heat is being added to the fluid at a point where its 
velocity equals the local velocity of sound (M «= 1), and that the 
maximum, temperature of the fluid occurs at M * l/v'“7" where 7 is 
the ratio of specific heats. These results are illustrated, by 
application to high-speed flow through combustion chambers. 


INTRODUCTION 

The study of combustion phenomena in a flowing gab must include 
consideration of the thermal and the dynamical properties of the 
fluid as well as of its chemical nature. The flow pattern of a 
burning gas, for example , is dependent upon the quantity of heat 
liberated at each inter nal point by the combustion of the moving 
fluid. 


The purpose of this report is to summarize, without extended 
proofs, the results of a study of a simplified model of nonadiabatio, 
compressible fluid flow, both subsonlo and supersonic, and to state 
these results in a form that will make them immediately useful in 
providing a theoretical background for current technical problems of 
high-speed combustion. The model chosen is the flow, in a tube of 
constant cross section, of a compressible fluid to which heat is 
being added. This model is a natural generalization of fluid-flow 
models commonly used and offers a logical point of departure for 
more advanced calculations. 

In the present treatment, the local Maoh number M of the 
flow is extensively used as a variable and the effects of heat addi- 
tion upon the limiting conditions of flow of a compressible fluid 
are emphasized. Apparently these effects have not previously been 
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explicitly described in the literature in terms of M, although 
related subjects in one-dimensional compressible flow have been 
treated using U by Bailey (reference 1) and Nielsen (reference 2) . 


THEORY 

Statement of the Basic Equations 

The characteristics of the model of nonadiabatic compressible 
flow, which were chosen in order to restrict and simplify the treat- 
ment, may be stated explicitly as follows: 

1. The flow is steady (independent. of time). 

2. Pressure, temperature , and velocity do not vary within any 
section normal to the flew. The rate of heat addition to the fluid 
is therefore constant within a section. 

3. The coefficients of viscosity and of thermal conductivity 
are equal to zero; hence, the fluid flow Is not affected by fric- 
tion or by conduction of heat in the direction of flow. 

U. The fluid obeys the thermal equation of state far a perfect 

gas. 

% . The values of the specific heats and of the gas constant 
are invariant throughout the f lev/ . 

With these restrictions, the fluid must move in accordance with 
the following equations, which describe its nonadiabatic fretion- 
less flow in a pipe (fig. l) : 


T 2 I 2 

Cp Ti + + Q - Cp T 2 + 7^2 (conservation of energy) 

n p 

t>i + Pi V - p 2 + P 2 Vg (equation of motion) 

- p£ Vg (conservation of mass) 


Pl T]_ p 2 T 2 


( 1 ) 

( 2 ) 

(3) 

U) 


(thermal equation of state) 
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The symbols p, p f T, V, and Cp In these equations represent 

pressure, density, temperature, velocity, and specific heat at 
constant pressure with consistent unite used- throughout . The sub- 
scripts 1 and 2 indicate that the symbols refer to sections 1 
and 2 of the flow pictured In figure 1. The symbol Q represents 
the heat added per unit mass of fluid between sections 1 and 2. 

A more convenient form of the equations is obtained by intro- 
ducing the local Mach number M « 7/y ^RT in place of the velocity 
where R 1 b the gas constant. Equations (1), (2), and (5) then 
become 

T i ( x ♦ q + “i 2 ) ■ T z f + “a 2 ) <5) 

Pl (l + Y Mi 2 ) - p 2 (l + Y Mg 2 ) (6) 

f>i a ^2 *t/r 2 (7) 

The quantity q is equal to Q/(CpT^ ) .'ind Y is oqual to Cp/Cy. 

Through multiplication and division of cqaations (5), (6), and 
(7) it is possible to derive the fundamental equation that describes 
the addition of heat to n compressible fluid in one-dimensional 
motion . 



Application of the Second Law of Thermodynamics 


An important characteristic of the model described by the basic 
equations can be derived directly from tiro Second Law of Thermo- 
dynamics. The change in entropy per unit mass AS of a perfect gas 
between the two thermodynamic states in which it exists at sec- 
tions 1 and 2 (fig. 1) is given by the expression 


AS 


1— Y 



(9) 


h 
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By use of equations (U), (6), and (7) the 
of the fluid in passing from section 1 to 
in terras of the Mach numbers and M 2 


specific-entropy change 
section 2 may be rewritten 
as 


r+1 


AS - 



( 10 ) ■ 



+ r k 


2 \ 

*2 J 


Y 


It is s cr-n that the entropy change of the moving fluid depends 
only upon the Mach number. Detailed examination of equation (10) 
discloses the fact that the entropy when considered as a function 
of M 2 possesses only one maximum, which occurs at = 1. 

’Then Up a 1, therefore, the entropy caimot increase. Because 
the entropy change of a system that behaves in accordance with the 
Second Law of Thorraudynamics must be positive when the quantity of 
heat added to it is positive, it is concluded that addition of a 
positive amount of heat to the fluid is impossible when M 2 ■ 1. 


RESULTS AND DISCUSSION 

From equation (8) and previous equations, a number of detailed 
characteristics of the steady nonadialr atic flow of a gas in a pipe 
of constant section (fig. 1) may hi derived. 


Mathematical Description of the Flow 


The heat added q » Q/(CpTq) and the changes in pressure, 
temperature, density, and velocity between sections 1 and 2 may- 
each be expressed as a simple rational function of M]_ and M 2 ,' 
as was done for q in equation (8). One physically significant 
value of M 2 exists for each and q > 0 ii M]_ < 1 and two 


such values if tlj > 1 . 


real value of r.lpl if 


If Ki< 


iS (Yq + 1)* 


Y - 1 


1 >Mj : 




Y - 1 


2Y (Yq + lj* 


there is only one 


there are two real 


values of V^, one of which corresponds to an expansion shock and 
may be shown to be impossible on thermodynamic grounds. For M^>1, 
transition from th-> larger of the two possible values of M 2 to tho 
smaller one corresponds to the formation of a compression shock. 
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The direction of change of the quantities p, p, T, and Q 
between sections 1 and 2 can be specified by the signs of the deriv- 
atives of . these quantities : with respect to M. For the entire range 
of M from zero to infinity 


& < 0 

dM - 


^>0 

dM “ 


l£< o 

dM ■ 


a - V 2 )' 1 £ > O 




(ii) 


As heat is added to the compressible fluid, the Mach number 
■will increase for all values of M <3j on the other hand, the tem- 
perature will rise only for M< l/\/Y or M >1 but will fall in 
the range of Mach number l>M>l/vY. This behavior can be derived 
directly from equations (11). 


limiting FIovt Conditions 


If the values of given by equation (8) are to be real, 

q must be limited and therefore the heat that can be added to a 
fluid moving with Mach number Mi is limited* the limit is 
expressed by the inequality J . 2 


< ( Bl ~ 

2 (T + 1) 


( 12 ) 


Thus, in agreement with the conclusion based on the Second Law of 

Q 

Thermodynamics, the heat added q =• 


*1 


cannot be positive for 


M x - 


1, that Is, for flow velocity equal to the local velocity of 
sound. The quantity of heat indicated by the equality sign in 
equation (12) is the quantity that is necessary to increase or 
decrease the Mach number to unity from its initial value Mi» If 
heat is added to but not removed from the fluid, this value of 
unity for the Mach number will occur at the downstream end of the 
tube, provided that the heat is supplied entirely by heat transfer 
from the wall of the tube. 


The limit to heat addition for a given entrance Mach number 
may also be thought of as a limit to the entrance Mach number for 
a given heat addition. The effect of heat addition upon the limit 
to entrance Mach number is represented In figure 2 by the curve 
for Mg » 1. 
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Because the mass flow Is proportional to (MipiJ/^/Ti, the 
limit to heat addition may be alternatively described as a limit 
to the mass flow for a given amount of heat added. If M^< 1 
and M does not exceed unity between sections 1 and 2, there will 
be an upper limit to the mass flow, as discussed in reference 3» 

It may be shown that the same upper limit exists for Mi < 1 if li 
does exceed unity between sections 1 and 2. If Mi>l, a lower 
limit to the mass rate of flow without shock exists, which is the 
same whether or not M is everywhere greater tlian unity between 
sections 1 and 2. The possibility of continuous transitions from 
1 to K >1 and vice versa are discussed in the next section. 

In all instances the limiting flow conditions can be most easily 
expressed analytically in terms of Mi as the independent varible. 
There is greater physical reality, however, in the expression of 
the limiting flow conditions in terms of a limit on the mass flow. 
Consequently Mi, which always bears a simple relation to the 
mass flow, has been taken to be the dependent variable in figure 2 
and in other similar figures. 

Tg (l + Y Ml 2 ') 2 

The temperature ratio — - cannot exceed 811(1 

T 1 hr Mi 

attains this maximum value at Mg <* l/VY* for all values of Mi 
and q. This limitation can be expressed alternatively as a limit 
to the entrance Mach number, which is then a function of Tg/Ti 
alone and which occurs when Mg ■ l/v^ • A graphical representa- 
tion of tne limiting Mi is given in figures 2 and 3 for both 
Mg - 1/y/T and Mg = 1. If the temperature ratio Tg/^l ^- s sp' 3C “ 
ified, the maximum entrance Mach inuribar (Mi < 1) is larger for 
Mg - l/VT” than for Mg =* 1. If q rather than Tg/Ti is spec- 
ified, the limiting entrance Mach number (Mi < 1) is less for 
Mg - l/v^ than for Mg - 1, as is evident from figure 2 -and 
the previous discussion. 

The minimum value of the pressure ratio Pg/Pi corresponding 
to the maximum allowable heat addition and to Mg ■ 1 is equal to 
(l + Y Mi) /(l + Y) for subsonic flow. If the pressure ratio is 
decreased below this limit by decreasing p ? , the flow will not 
change because the effects of the pressure f 'eduction cannot be 
propagated upstream through the region where Mg = 1 nt the down- 
stream section 2. tor supersonic flow, the same expression gives 
the maximum pressure ratio that can occur with steady flow. 
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Transition between Subsonic and Supersonic Flow 

The theory indicates that a continuous transition from sub- 
sonic to supersonic' flow or vice versa can be effected by addition 
of heat until !t « 1 followed by removal of heat. If heat iB 
added to but not removed from the moving fluid, such a continuous 
transition is not possible. A discontinuous transition from sub- 
sonic to supersonic flow (adiabatic or nonadlabatic expansion 
shock) may be shown to be impossible on the basis of the Second 
Law of Thermodynamics. A discontinuous transition from supersonic 
to subsonic flow is, however, always possible through the action 
of a compression shock, whether the heat added is positive, zero, 
or negative. The action of viscosity or of thermal conductivity, 
which are omitted from consideration in the present treatment, 
would provide the mechanism whereby such a discontinuous transition 
would be effected (reference li). 


Thermal Choking 

It has been shown that steady flow cannot occur unless the flow 
conforms to limitations on the pressure and temperature ratios. 

Iffhen conditions that exceed these limitations are initially imposed 
on the flow, the flow will adjust itself in such a way as to again 
reach conformity with the limitations, 'lhis self -adjustment of the 
Clor is somewhat similar to the adjustment that occurs when a per- 
fect compressible fluid flows (1) isentropically through a tube of 
varying section, (2) adiabatically with frictional pressure drop 
through a tube of constant section (reference or (3) iso ther- 
mally with frictional pressure drop through a tabu of constant 
section (reference 6). By analogy with these types of flow, the 
term "thermal choking" might be appropriately applied to the phe- 
nomenon of limited mass flow discussed in the present paper. The 
similarity between the various types of choking phenomenon is being 
investigated in more detail. 


APPLICATIONS 

The theoretical results indicate that the mass flow of a com- 
pressible fluid cannot exceed a limit that Is fixed by the values 
of p-^, T^, and q. For each value of Mj or for given values 

of p-^ and T^, at the first section, there is only one value 
of q or of Q, which corresponds to the limiting flow condition. 
These results can be applied to the flow of a burning mixture 
through the combustion chamber of a Jet-propulsion power plant. 
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It is assumed that the cross sections of the ducts ■which 77111 
be considered are constant and that frictional effects, change of 
gas constant and of Y can be neglected. The maximum mass flow C 
for subsonic velocities throughout the chamber has been computed 
for the inlet conditions ■ 15 pounds per square inch and 


500° F absolute and is plotted in figures ii, 5» and 6 for 


M 2 “ l/y^Y and M 2 “ 1. The curves in figures U 


s, 

and 


5 were 


obtained from figures 2 and 5, respectively, by changing the scales 
in accordance with the chosen entrance conditions. The data in 

Mj 2 
V 


figure 6 were derived from the expression 


unions, ine 

\pj- 


which 


is discussed in a previous section. Because the fuel-air ratio 
determines the value of Q, the heat released per pound of mixture 
for complete combustion of a given fuel#points can be specified on 
the curves of maximum mass flow that correspond to a range of values 
of fuel-air ratio. Although the very lean mixtures corresponding 
to fuel-air ratios less than about 0.02 would not. burn under ordi- 
nary circumstances, the values of fuel-air ratio can still be used 
as a convenient measure of heat release. In figurer h, 5, and 6 
the fuel-air ratios are indicated for a fuel whose heat of combus- 
tion is 18,000 Btu per pound. 


The determination of the conditions for maximum flow for com- 
plete combustion (Mg - 1 ) will servo to illustrate the use of 
the figures. If the limit on the total gas outflow (including 
secondary air) from the combustion chamber is desired, an over-all 
fuel-air ratio of about 0.02 is appropriate. The limit on the gas 
flow through a burner duct containing a richer mixture with fuel- 
air ratio equal to 0.05 is calculated for comparison. From the 
curve for ilg ■ 1 (fig. U), the maximum flow rates for these two 
Cuel-air ratios are found to be 22.2 and lit . 8 pounds per square 
foot per second, respectively. The final temperatures determined 
from figure 5 are 1630 ° and 3 WiO° Tr absolute, and the final pres- 
sures from figure 6 are 6.80 and 6 ,ii 8 pounds per square inch. The 

Q 

entrance Mach numbers » ftc >~0 corresponding to these flow rates 


are 0.250 and 0.166 or, because the velocity of sound at the entrance 
temperature of 500° F absolute is 1096 feet per second, the entrance 
velocities arc 27h and 182 feet per second. The values of can 

Q Q 

be chocked by calculating the values of q « ^ =■ =- 5 -?, which are 

Cp ■. It! U 

2.9U and 7.1U, and by reading the values of Mj (fig* 2). It is 
evident that the mass flow of burning gasos can be considerably 
restricted by the effects of thermal choking. 
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If tbs upstream pressure Is Increased, above 15 pounds per 
square inch, the -downstream pressure and the flow should Increase 
proport i 05 sately with the pressure, because Mj_ and Mg remain con- 
stant, If the downstream pressure Is decreased while remains 
equal to 15 pounds per square Inch, then Mj_, Mg, and the flow 
should each remain cons tent. The theoretical limits placed by the 
occurrence of thermal choking an the flow are therefore similar 
to those found for Isentroplo flow from converging nozzles. 


Aircraft Engine Resear oh Laboratory, 

Rational Advisory Committee for Aeronautics, 
Cleveland, Ohio. 
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Fig. I 


Pl» Pl» T 1 


P2» P2* t 2 




Q=heat added per 
sections 1 and 


unit mass between 
2 
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Figure 1. - One- dimensional fluid flow with heat 
addition. 








Entrance Mach number 


Fig. 3 
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Maximum mass flow, <J, lb/(sq ft) (see) 



NACA ACR No. ESA29 


Maximum mass flow, G, lb/(sq ft) (sec) 





Exit temperature, Tq> °P absolute 

Figure 5. - Relation between maximum mass flow and exit temperature for 
entrance pressure of 15 pounds per square inch and entrance temperature 
of 5000 p absolute. 
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Maximum mass flow* G* lb/(sq ft} (sec) 
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Fig. 6 
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exit pressure for entrance pressure of 15 pounds 
per square inch and entrance temperature of 
500° F absolute. 




